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Abstract 

The framework of multifractal analysis (MFA) is distilled to the most sophisticated 
one. Within this transparent framework, it is shown that the harmonious represen- 
tation of MFA utilizing two distinct Tsallis distribution functions, one for the tail 
part of probability density function (PDF) and the other for its center part, explains 
the recently observed PDF's of turbulence in the highest accuracy superior to the 
analyses based on other models such as the log-normal model and the p model. 
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1 Introduction 



The multifractal analysis (MFA) is a unified self-consistent approach dealing 
with those systems exhibiting large deviations. It has been proposed by the 
present authors [1,2,3,4,5,6,7,8,9] in order to analyze various probability den- 
sity functions (PDF's) observed in the system of fully developed turbulence. 
The PDF's, given for variables normalized by their own standard deviations, 
represent a characteristic large wings (or fat tails) due to the intermittency 
specific to turbulence. Since the pioneering work by Kolmogorov [10], there 
appeared so many attractive attempts to analyze the large wings, e.g., the log- 
normal model [11,12,13], the p model [14,15] etc., but no one was satisfactory. 
The MFA is constructed with the help of the Tsallis-type distribution func- 
tion [16] that provides an extremum of the extensive Renyi entropy [17] or the 
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non- extensive Tsallis entropy [16,18] under appropriate constraints [3]. After 
a rather preliminary investigation of the p model [1], we developed further to 
derive the analytical expression for the scaling exponents of velocity structure 
function [2,3], and to determine the PDF's of the velocity fluctuations [3,4,5], 
of the velocity derivative [6] and of the fluid particle accelerations [7]. 

In this paper, we will show that a harmonious representation of PDF's by 
means of two distinct Tsallis-type distributions provides us with a descrip- 
tion for experimentally or simulationally observed PDF's in the highest accu- 
racy compared with other multifractal model such as the log-normal model 
[11,12,13] and the p model [14,15]. In order to perform a transparent compar- 
ison, the framework of MFA is rephrased in its most sophisticated fashion in 
which the tail part of PDF, giving the probabilities of events larger than its 
standard deviation, is written down once the multifractal spectrum for the spa- 
tial distribution of singularities is specified (see (1) below), whereas the center 
part, giving the probabilities of events smaller than its standard deviation, is 
assumed to be analyzed by the Tsallis-type PDF for the variable itself [8,9] 
(see (2) below). The various PDF's extracted out from the DNS conducted by 
Gotoh et al. [19,20], and the PDF of fluid particle accelerations observed by 
Bodenschatz et al. in their Lagrangian measurement of particle accelerations 
[21,22,23] will be analyzed with the help of the multifractal spectrums for the 
log-normal, the p model and the harmonious representation. 

MFA of turbulence rests on the scale invariance of the Navier-Stokes equation 
for high Reynolds number, and on the assumption that the singularities due to 
the invariance distribute themselves, mult ifr act ally, in physical space. Let us 
consider the fluctuation 5x n = \x(»+£ n )—x(») \ of a physical quantity x at the 
nth multifractal step satisfying the scaling law \x n \ = \5x n /5xo\ = 5 n °"^^ 3 with 
$n — £ n /£o = &~ n (n — 0, 1, 2, • • •). We call n the multifractal depth which can 
be real number in the analysis of experimental data. In the following in this 
paper, we will put 6 = 2 that is consistent with the energy cascade model. At 
each step of the cascade, say at the nth step, eddies break up into two pieces 
producing the energy cascade with the energy-transfer rate e n that represents 
the rate of transfer of energy per unit mass from eddies with diameter £ n 
to those with £ n+ i. Then, we see that the derivative \x'\ = lim^^x^ with 
the nth difference x' n = 5x n /£ n for the characteristic length £ n diverges for 
a < 3/0 [24]. The velocity fluctuation 8u n is given for = 1, and the pressure 
fluctuation 5p n is for = 2. The velocity derivative and the fluid particle 
acceleration are defined, respectively, by \u'\ = lim^oo u' n and |a| = lim^oo a n 
with the nth velocity difference u' n = Su n /£ n and the nth pressure difference 
(the nth acceleration) a n = 5p n /£ n . Note that the fluid particle acceleration 
a is given by a = du/dt + (u ■ V)n. For the energy transfer rate, = 3 since 
it has the scaling relation e„/e = Sj? -1 - 

The real quantity a is introduced in the scale transformation [25,26] x — > x' — 
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Xx, u -> v! = \ a / 3 u, t -> f = A 1_a/3 t, p -> pt = \ 2a / 3 p that leaves the 
Navier-Stokes equation du/dt + (w • V)w = — Vp + z/V 2 w of incompressible 
fluid invariant when the Reynolds number Re = <5wi n £i n /z/ is large. Here, v is 
the kinematic viscosity, and p = p/p with the thermodynamical pressure p and 
the mass density p. 5u in and t in represent, respectively, the rotating velocity 
and the diameter of the largest eddies in turbulence. The largest size of eddies 
is, for example, about the order of mesh size of a grid, inserted in a laminar 
flow, that produces turbulence downstream. 



2 General framework 



MFA starts with an assignment of the probability, to find a singularity speci- 
fied by the strength a within the range a ~ a + da, in the form P( n \a)da = 
(\f"(oto)\ | ln<5 n |/27r) 1 / 2 6^~^°^ da [15,3]. Here, f(a) represents an appropriate 
multifractal spectrum defined in the range a m i n < a < a max - Note that f(a) 
does not dependent on n because of the scale invariance. The multifractal spec- 
trum is related to the mass exponent r(q), defined by ((e ra /e) 5 ) = a^g <5~ T ® +1 ~^ 
with a>3q = (|/"(ao)|/|/"( Q; g)l) 1 ' /2 ) through the Legendre transformation [15]: 
f(a) = aq + r(q) with a = aq = —dr(q)/dq and q = df(a)/da. The average 
(• • •) is taken with P( n \a), and a = aq =0 = (a). 

The PDF n^(x n ), normalized as / dx n Tl^\x n ) — 1, is assumed to consists 
of two parts, i.e., U.^\x n ) = U^(x n ) + AH < ^'\x n ) where the first term is 
related to P^(a) through U { ^(\x n \)d\x n \ = [(1 - 27$) /2}P^\a)da with 
the transformation of the variables \x n \ = 5 n a ^^ 3 , and the second term is re- 
sponsible to the contributions coming from the dissipative term in the Navier- 
Stokes equation violating the invariance under the scale transformation given 
above. Here, we introduced 7^ by 27^ = / dx n \x n \ m AU ( ^'\x n ). Note that 
we are dealing with the symmetrized part of PDF's by assuming that the 
large deviation stemmed from the singular first term contributes to symmet- 
ric part of PDF's. We have the rath order structure function for x n in the 
form i\x n \ m ))4 = J dx n \x n \ m uf ] (x n ) = 2 7 g, + (1 - 2 1 $)a (f)m 5n 4 ' rn with the 
corresponding scaling exponents ( 3q - = 1 — r(q). For 0=1, C<f>m reduces to the 
scaling exponents ( m of the mth order velocity structure function. 

Let us introduce the PDF f[^(£ n ), to be compared with observed data, defined 
by n^ n) (£ n )d£ n = Tl { £\x n )dx n or by n^ n) (£ n )d£ n = Yl ( f' \x' n )dx' n for the vari- 
able £ n = x n / ((xl)) 1 / 2 = x' n / {{(x'n) 2 )) 1 ^ 2 both of the fluctuation x n and of the 
derivative x' n normalized by their own standard deviations. For values of the 
variable |£ n | larger than the order of its standard deviation, £* < |£ n | < £™ ax 
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(equivalently, a min < |a| < a*), the PDF is given by [8,9] 

n^ia^nl = U^(\x n \)d\x n \ = [(1 - 27$)/2]P (B) («)<k (1) 



where £•« = ^/ Q - n/3 -^ /2 with & = [2 7 g^ + (1 - 2 7 $)o 2 ,]- 1 /*. 
Note that £* ~ 1 as can be seen below when we analyze experiments. This 
tail part represents the large deviations, and manifests itself the multifractal 
distribution of the singularities due to the scale invariance of the Navier-Stokes 
equation when its dissipative term can be neglected. The entropy index q 
should be unique once a turbulent system with a certain Reynolds number is 
specified. For smaller values of the variable, |£ n | < £* (equivalently, a* < \a\), 
we assume that the PDF has the Tsallis-type structure with a new entropy 
index q' [8,9] 



n£°(£nK„ = lnSS(s B ) + Anjr'(s B )| dx n 



r(«) 



n 



(n) 



{i - a - q ') (0 + 3/v)) [(ucf - 1] /20} 1/(1 ^' ) de„ 



(2) 



where IT^ = 3(1 - 2 7 g)|/' / (a o )| 1/2 /20a(27r| lnt^) 1 / 2 . This center part is 
responsible to smaller fluctuations of the variable, compared with its standard 
deviation, stemmed from the dissipative term violating the scale invariance. 
The entropy index q' can be dependent on the distance of two measuring 
points. 

The two parts of the PDF, (1) and (2), are connected at f * = £ n <5/ Q * /3_C ^ /2 
under the conditions that they have the common value, n^ n \ and the common 
log-slope, —(0 + 3/'(a:*))/0£* there. The value a* is the smaller solution of 
C20/2 — 0a/3 + 1 — f(a) = 0. The point £* has the characteristics that the 
dependence of n^(£*) on n is minimum for large n (see Fig. 4 in [4]). With 
the help of (1) and (2), we obtain AII^(x n ), and have the analytical formula 



27$, = TO - £JK,)/(1 + - where 



■(n) 



(n) 



<l = ( 3 e +1)0aV3 -^ /2 /0) Vl/"WI/27r|ln(y n | 

l-(l-g') (0 + 3/V))(^ 2 -l)/20; 

^max 

4^ = ^V / l/ / '(«o)||ln^|/2 7 r y da^/»-/W. 



1 

x y ^ [z[ m 




1/(1-9') 



(3) 
(4) 



For later convenience, we introduce here the quantity £ n = £ n <5n' 



(pao/3— C2</>/2 
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The PDF's rT n ) (£ n ) of the velocity fluctuations and of the velocity derivatives 
are give by the common formula & n \C, n ) = fl^ 1 (^ n ) for the normalized vari- 
able £ n = 5u n /{{(8u n ) 2 )) 1 ' 2 , whereas the PDF's k (n \u n ) of the pressure fluctu- 
ations and of the fluid particle accelerations are given by the common formula 
A.( n \u n ) = n^ 2 (o; n ) for the normalized variable uo n = 5p n / (((Spn) 2 }} 1 ^ 2 . 



3 Log-normal model 

In the log-normal model [11,12,13], one consider the ratio e n je n ~\ (n = 1, 2, • • •) 
as independent stochastic variables, and apply for n ^> 1 the central limit 
theorem to the summation of their logarithms, (W 2 )" 1 / 2 Z)™=i m ( e j/ e j-i) — 
(no- 2 ) -1 / 2 ln(e n /e) = ^Jna^ 1 ^— a) In 5, to have the Gaussian distribution func- 
tion P {n \a) = (n/27ta 2 ) 1 / 2 e- nia - ao)2/2cT2 for the range -oo < a < oo. Here, 
we used the scaling relation between e n and a. Then, we have the multifractal 
spectrum and the mass exponent in the forms f{a) = 1 — (a — a§) 2 j2<j 2 In S 
and r(q) = 1 — a q + g 2 <r 2 (In 5) /2, respectively. We see that a 3q - = 1. 

The dependence of the parameters «o an d ° on the intermittency exponent 
/i is determined with the help of the two independent equations, i.e., the 
energy conservation: (e„/e) = 1 (equivalently, r(l) = 0), and the definition 
of the intermittency exponent /i: (e 2 /e 2 ) = <5~ M (equivalently, /i = 1 +r(2)). 
Here, e is the energy input rate to the largest eddies. The parameters are 
specified by means of fi as a = 1 + /i/2 and a 2 = /i/ln<5. Then, we have 
f{ot) = 1 — (a — ao) 2 /2fi and r(q) = (1 — q)Dq with the generalized dimension 
Dq = 1 — /ug/2, which are the same as derived in [15]. We know that aq = 
a - fiq. 

The explicit form of the PDF (1) for £* < |£ n | < oo is found to be 

(5) 



n Un) = n^— exp 



(31n|a/£n,o|) 2 



20 2 /i| ln5 n | 

The connection point is given by a* = ao — (V3 — l)0/x/3. 
4 P model 



The distribution function P^ n \a) for the p model [14,15] is specified based on 
the binomial multiplicative process in the form [1] P {n \a) = (4 n) )- 1 [2^(l - 
y) 1 ~ v \~ n with y = y(a) — [a + log 2 (l — p)]/log 2 [(l — p)/p\ and the partition 
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function Zq = y 7r/2nlog 2 [(l —p)/p] for n»l. The multifractal spectrum 
is given by f(a) = -{y(a) log 2 y(a) + [1 -y(a)] log 2 [l -y(a)]}, which leads to 
the mass exponent r(q) = log 2 [p q + (l— p) q }. We have f"{otq) = — ln2/{ln[(l — 
p)/p}} 2 yq( 1 - Vq) wi th yq = y(a q ) = p q /[p q + (1 -p) 9 "], a q = -[p q ~\og 2 p + (1 - 
p) 9 "log 2 (l -p)}/{p q + (1 -p) 9 "], and then we obtain a 35 = 2[p«(l - pf} 1 ' 2 /[p^ + 
(1 - pf}. We see that, for p > 1/2, a min = - log 2 p, a max = - log 2 (l - p). 

The dependence of p on the intermittency coefficient \i is derived through its 
definition /i = 1 + r(2) to give p = (1 + - l)/2. 

The explicit form of the PDF (1) for £ < |f n | < £™ ax reduces to 



= 1 + 2 (f - ?/o)] 1+2( ^ o) [1 - 2 (y - y,)] 1 - 2 ^}^ (6) 

lsn| 

with y - y = 31n(|a|/|a,o|)/0{log 2 [(l - p)/p}}\ \n6 n \. 



5 Harmonious representation 

In the harmonious representation of MFA, we adopt the Tsallis-type dis- 
tribution function [2,3] P^ n \a) = (zf )" 1 -^ - [(a - a ) /Aa] 2 }^^ with 
(Aa) 2 = 2X/[(l—q) ln2] and the appropriate partition function Z^K Here, q is 
the entropy index introduced in the definitions of the Renyi [17] and the Tsallis 
entropies [16,18]. This distribution function provides us with the multifractal 
spectrum f(a) = 1 + (1 — q)^ 1 log 2 [l — (a — ao) 2 /(Aa) 2 ] which, then, produces 
the mass exponent r(q) = 1 — a^q + 2Xg 2 (l + C^ 2 ) -1 + (1 — — log 2 (l + 

C l q ' 2 )} with C q =l + 2q 2 {l - q)X In 2. We have f"(a q ) = -C l J 2 (\ + C 1 q /2 )/2X 
with a- q = a - 2qX/(l + C 9 - /2 ), hence a 3q = {2/[C l q /2 (l + C- 72 )]} 1 / 2 . We see 
that a min /max = a T Aa. 

The dependence of the parameters ao, X and q on the intermittency expo- 
nent fj, is determined, self-consistently, with the help of the three independent 
equations, i.e., the energy conservation: r(l) = 0, the definition of the inter- 
mittency exponent \i: /i = 1 + r(2), and the scaling relation: [30] 1/(1 — q) = 
l/a_ — l/ct + with a± satisfying f(ct±) = 0. Within the range 0.13 < // < 0.40 
where most of the experiments are covered, the self-consistent equations are 
solved, numerically, to give [6] ao = 0.9989 + 0.5814/i, X = 1.198// and 
q = -1.507 + 20.58/i- 97. 11/i 2 + 260.4//- 365.4// + 208.3/A It is interesting 
to compare the first two with the corresponding equations for the log-normal 
model (q — 1), i.e., a = 1 + 0.5// and X — //. 
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Fig. 1. Analyses of the PDF's of the velocity fluctuations (closed circles) for three 
different measuring distances, observed by Gotoh et al. at R\ = 380, with the 
help of the PDF's ]T n )(£ n ) by the harmonious representation (solid line) and by 
the log-normal model (dashed line) are plotted on (a) log and (b) linear scales. The 
PDF's by the p model (dotted line) are compared with the PDF's by the harmonious 
representation (solid line). Comparisons are displayed in pairs. The solid lines in 
each set of pairs are the same. For better visibility, each PDF is shifted by —2 unit 
in (a) and by —0.2 in (b) along the vertical axis. Parameters are given in the text. 

The explicit form of the PDF (1) for f * < |f n | < £™ ax turns out to be 



ni%n) = ii7^ 



r (n) £,n 



1 - 



q (3]n|£ n /f n>0 |)' 



n 



20 2 X|ln5 n | 



n/(l-q) 



(7) 



6 Analysis of experiments 



The scaling exponents ( m of velocity structure function reported by Gotoh et 
al. [19] are analyzed by the method of the least squares (MLS) with the theo- 
retical formulae of the harmonious representation, of the log-normal model and 
of the p model, giving, respectively, the values of the intermittency exponent 
fi = 0.240, 0.217 and 0.249. 

In Fig. 1, the PDF's of the velocity fluctuations (closed circles) measured by 
Gotoh et al. in their DNS at R\ = 380 [19] for three different measuring 
distances, r/r] = £ n /i] = 2.38, 19.0, 1220 from the top set of pairs to the 
bottom set, are analyzed with the help of the PDF's IT") (£ n ) of the harmonious 
representation (solid line) and of the log-normal model (dashed line). Here, 
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Fig. 2. Analyses of the PDF's of fluid particle accelerations, measured by Gotoh 
et al. at R\ = 380 (circles in the top set) and by Bodenschatz et al. at R\ = 690 
(circled in the bottom set), by means of the PDF's M n \uj n ) by the harmonious 
representation (solid line) and by the log-normal model (dashed line) are plotted on 
(a) log and (b) linear scales. The PDF's by the p model (dotted line) are compared 
with the PDF's by the harmonious representation (solid line). Results are displayed 
in pairs. The solid lines in each set of pairs are the same. For better visibility, 
each PDF is shifted by —2 unit in (a) and by —0.4 in (b) along the vertical axis. 
Parameters are given in the text. 

7] = (z/ 3 /e) 1//4 represents the Kolmogorov scale. The DNS data points are 
symmetrized by taking averages of the left and the right hand sides data. In 
each set, a comparison of the PDF of the harmonious representation (solid 
line) with the PDF of the p model (dotted line) is given. For the harmonious 
PDF's (solid line), q = 0.391 (p = 0.240), and, from the top set to the bottom 
set, (n, q') = (20.7, 1.60), (13.6, 1.50), (6.10, 1.20), Q = 1.10, 1.23, 1.43 {a* = 
1.07) and C ax = 204, 38.2, 6.63. For the PDF by the log-normal model (dashed 
line), from the top set to the bottom set, (n, q') = (21.5, 1.70), (13.0, 1.63), 
(5.00, 1.24) and = 1.19,1.34,1.51 (a* = 1.06). For the PDF by the p 
model (dashed line), from the top set to the bottom set, (n, q') = (21.0, 
1.60), (13.0, 1.62), (5.50, 1.20), £* = 1.06,1.07,1.55 (a* = 1.08) and £™ x = 
19.3, 7.70, 3.31. The tail of PDF of the p model stops at £™ ax whi ch is smaller 
than the maximum value of measured data point for each measuring distance. 
Note that every PDF's are plotted on (a) log and (b) linear scales. 

In Fig. 2, the PDF's of the fluid particle accelerations (closed circles) reported 
by Gotoh et al. at R\ = 380 (top set) [20] and by Bodenschatz et al. at 
R x = 690 (bottom set) [21,22,23] are analyzed with the PDF's A (n) (cu„) of the 
harmonious representation (solid line) and of the log-normal model (dashed 
line). The measured data points both on the left and right hand sides of the 
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PDF's are shown altogether on one side by closed circles in the figure. In 
each set, a comparison of the PDF of the harmonious representation (solid 
line) with the PDF of the p model (dotted line) is given. For the harmonious 
PDF's (solid line), from the top set of pairs to the bottom set, q =0.391, 0.391 
(pi = 0.240, 0.240), (n, q') =(17.5, 1.70), (17.1, 1.45) w* = 0.622, 0.605 (a* = 
1.01, 1.01), and < iax = 2530, 2040. For the PDF by the log-normal model 
(dashed line), from the top set to the bottom set, (n, q') = (17.0, 1.50), (18.5, 
1.04), u* = 0.644, 0.558 (a* = 1.00, 1.00), and u™ x = 49.6, 76.6. For the PDF 
by the p model (dotted line), from the top set to the bottom set, (n, q') = 
(19.0, 1.50), (18.5, 1.20) u* = 0.539, 0.547 (a* = 1.01, 1.01), and u™ x = 54.1, 
45.2. Note that every PDF's are plotted on (a) log and (b) linear scales. 

The values of n for all the PDF's are determined by MLS by adjusting the 
integrand £^n^(£„) of the fourth moment both of data and of the theories. 
The values of q' for all the PDF's are obtained by MLS by fitting the theoretical 
PDF's at the center part, |£„| < £*, with the observed PDF's. 



7 Discussions 



As was shown in this paper, the harmonious representation within MFA pro- 
vides the highest accuracy in analyzing observed PDF's. This high accuracy 
allows us to extract some information about the dynamics underlying turbu- 
lence. For example, through the analysis of the velocity fluctuations in Fig. 1, 
we extracted the dependence of q 1 on r/r) as [9] q' = 1.71 — 0.05 log 2 (r/?y). 
We then obtain the diffusion coefficient for the fluctuations of magnitude 
|£„| (< C) in the form D^ n )/ D^O) = 1 + [0.71 - 0.05 log 2 (r/r ? )][(0 + 
3/'(a:*))/20](£ n /£*) 2 . Note that this is dependent on the distance r of two 
measuring points. In deriving this, we assumed that the drift term is lin- 
ear since the magnitude of the fluctuations in the center part of the PDF is 
smaller than the order of their standard deviations (|£ n | < 1), and, therefore, 
used the formula for anomalous diffusion given in [29]. The investigation to 
this direction is in progress, and will be reported elsewhere. 
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